Role of initial data in spherical collapse 
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We bring out here the role of initial data in causing the black hole and naked singularity phases as 
the final end state of a continual gravitational collapse. The collapse of a type I general matter field 
is considered, and it is shown that given the distribution of density and pressure profiles at the initial 
surface from which the collapse evolves, there is a freedom in choosing rest of the free functions and 
dynamical collapse evolutions from this initial data, such as the velocities of the collapsing shells, 
so that the end state could be either a black hole or a naked singularity depending on this choice. 
Our results give some physical insights into these collapse outcomes as to how these phases come 
about. 
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An intriguing finding that has emerged from recent 
investigations on gravitational collapse scenarios in gen- 
eral relativity is that the final end state of such a collapse 
could be either a black hole (BH), or a naked singular- 
ity (NS). The theoretical and observational properties of 
these objects could be quite different from each other and 
so it is of immense interest to get an insight into how each 
of these phases come about as end states of a dynamically 
developing collapse ruled by gravity. When a star is suf- 
ficiently massive, and starts collapsing gravitationally on 
exhausting its nuclear fuel, it would not settle to a stable 
configuration such as a neutron star. What happens in 
such a case is an endless gravitational collapse ensues, 
where the sole governing force is gravity. Then either an 
event horizon of gravity forms well in advance to create 
a black hole, or in the case otherwise, the extreme den- 
sity regions may fail to be covered by the horizon, and 
a visible naked singularity may develop (see e.g. [1] for 
some recent reviews). Most of the models studied so far 
are spherical, and the matter cloud collapses under rea- 
sonable physical conditions such as an energy condition, 
and regularity of the initial data from which the collapse 
develops. The possible physical consequences of the later 
scenario has drawn some attention recently [2] . 

The major unresolved problem, however, is to under- 
stand what causes these BH/NS phases as end states of 
a continual collapse, and how these come about. In other 
words, we need to know why and when a naked singu- 
larity results as collapse end state as opposed to a black 
hole. We study here collapse for general Type I matter 
fields, and our purpose is to show explicitly that there ex- 
ist classes of solutions of Einstein equations which take 
the collapse evolution to either of the BH/NS phases, 
subject to usual regularity and energy conditions, and 
depending on the choice of the free functions available. 
While we discuss spherically symmetric collapse, the form 
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of matter we consider is rather generic which is any type 
I general matter field subject to an energy condition [3]. 
What is seen here is given the initial matter distribu- 
tion from which the collapse develops, there is a freedom 
to choose the velocity function (and its derivatives) for 
the collapsing shells, so that the collapse could end in ei- 
ther a black hole or naked singularity, depending on this 
choice. We thus get here some insight into what causes 
these phases, in terms of physical agencies such as ini- 
tial matter and velocity distributions, and the allowed 
possible dynamical collapse evolutions. 

The method we follow here is outlined below. In the 
case of a black hole developing as collapse end state, 
the spacetime singularity is necessarily hidden behind 
the event horizon of gravity, whereas in the case of a 
naked singularity developing, there are families of future 
directed non-spacelike trajectories coming out from the 
singularity, which can in principle communicate informa- 
tion to faraway observers in the spacetime. The existence 
of such families confirms the NS formation, as opposed 
to a BH end state. We shall study the singularity curve 
produced by the collapsing matter, the tangent to which 
at the central singularity at r — is related to the radi- 
ally outgoing null geodesies from the singularity if there 
are any, as we show here. By determining the nature of 
the singularity curve and its relation to the initial data 
and the classes of collapse evolutions, we are able to de- 
duce whether the trapped surface formation in collapse 
takes place before or after the singularity. It is this causal 
structure of the trapped region that determines the pos- 
sible emergence or otherwise of non-spacelike paths from 
the singularity. This settles the final outcome of collapse 
in terms of either a BH or NS. We point out that several 
familiar equations of state for which extensive collapse 
studies have been made so far, such as dust, matter with 
only tangential or radial pressures, and others form spe- 
cial cases of the considerations here. 

The spacetime geometry within the spherically sym- 
metric collapsing cloud can be described by the general 
metric in the comoving coordinates (t, r, 0, ef>) as given by, 
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ds 2 = -e 2 ^dt 2 + e 2 ^dr 2 + R 2 (t, r)dtf 



(1) 



where d£l 2 is the line element on a two-sphere. The 
energy-momentum tensor is that for any matter field 
which is type I (this is a broad class which includes 
most of the physically reasonable matter fields, includ- 
ing dust, perfect fluids, massless scalar fields and such 
others, see e.g. [3]), and is given in a diagonal form, 
T\ = -p,T 7 r = p r ,Tg = = p e . The quantities p, 
p r and pg are the energy density, radial and tangential 
pressures respectively. We take the matter field to sat- 
isfy the weak energy condition, i.e. the energy density 
measured by any local observer be non-negative, and so 
for any timelike vector V % , we must have TikV 1 V k > 0, 
which amounts to p > 0; p + p r > 0; p + pg > 0. 

Now for the metric Q the Einstein equations take the 
form, in the units (8irG = c = 1) 
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where we have defined G(t,r) 
e- 2v {R) 2 . 

The arbitrary function F = F(t, r) has an interpreta- 
tion of the mass function for the cloud, giving the total 
mass in a shell of comoving radius r. Energy conditions 
imply F > 0. In order to preserve the regularity at the 
initial epoch, we have F(ti, 0) = 0, that is the mass func- 
tion should vanish at the center of the cloud. As seen 
from equation (2), there is a density singularity in the 
spacetime at R = 0, and at R' = 0. However, the later 
one is due to shell-crossings and can be possibly removed 
from the spacetime [4], so we shall consider here only 
the shell- focusing singularity at R — 0, which is a phys- 
ical singularity where all matter shells collapse to a zero 
physical radius. We can use the scaling freedom avail- 
able for the radial co-ordinate r to write R = r at the 
initial epoch t = ti. Introducing the function v(t,r) by 
the relation, 



R(t, r) — rv(t, r) 



(6) 



we have v(U, r) 
0. The time t ■■ 



- l;v(t s (r),r) — and for collapse, v < 
t s (r) corresponds to the shell- focusing 



singularity R — 0, where all the matter shells collapse to 
a vanishing physical radius. 

From the point of view of dynamic evolution of the 
initial data prescribed at the initial epoch t = ti, there 
are six arbitrary functions of the comoving shell-radius r 
as given by, v(U,r) = v (r),ip(ti,r) = ip (r), R{U, r) = 
rp(U,r) = p Q {r),p r (ti,r) = p ro (r)p e (ti,r) = pg . We 
note that not all the initial data above are mutually in- 
dependent, because from equation @ we get, 
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Let us now assume that the initial pressures have phys- 
ically reasonable behavior at the center (r = 0) in that 
the pressure gradients vanish, i.e. p' r {0) — p'g o (0) = 0, 
and also that the difference between radial and tangential 
pressures vanishes at the center, i.e. p ro (0) — pg o (0) = 0, 
to ensure the regularity of the initial data at the center 
of the cloud. Then from equation©, it is evident that 
voir) has the form, 



u (r) = r 2 g{r) 



(8) 



where, g(r) is at least a C 1 function of r for r = 0, and 
at least a C 2 function for r > 0. 

Hence we see that there are five total field equations 
with seven unknowns, p, p r , pg, ip, v, R, and F, giving us 
the freedom of choice of two free functions. Their selec- 
tion, subject to the weak energy condition and the given 
initial data for collapse at the initial surface, determines 
the matter distribution and metric of the space-time, and 
thus leads to a particular time evolution of the initial 
data. 

Consider now a general mass function F(t, r) for the 
collapsing cloud, which can be written as 



F{t,r) = r 3 M{r,v) 



(9) 



Here M. is a regular and suitably differentiable function 
and from equations below, the regularity and finitcness 
of the density profile at the initial epoch t — ti requires 
that F goes as r 3 close to the center. The equations 
and (6) give, 
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Thus regular density distribution at the initial epoch is 
given by, po( r ) = 3M(r, 1) + rM(r, l) <r . It is evident 
that, in general, as v — > 0, p — > oo and p r — * oo. That 
is, both the density and radial pressure blow up at the 
singularity. 

What we show below is, given any regular initial den- 
sity and pressure profiles for the matter cloud from which 
the collapse develops, there always exist velocity profiles 
for collapsing matter shells, and classes of dynamical evo- 
lutions as determined by the Einstein equations, so that 
the end state of the collapse would be either a naked sin- 
gularity or a black hole, depending on the choice made. 
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In other words, given the matter initial data at the ini- 
tial surface t — t i: these evolutions take the collapse to 
end up either as a black hole or the naked singularity, 
depending on the choice of the class. 

To see this, we provide now classes of solutions to Ein- 
stein equations to this effect. Consider the class of veloc- 
ity profiles as determined (using equations (7) and (8)) 
by the general function, 



u(t,r) = A(t, R) 



(11) 



where A(t, R) is any arbitrary, suitably differentiable 
function of t and R, with the initial constraint A(ti, R) = 
vo(r). Using eauation (|ll|l in equation 0), we have, 
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(12) 



Here b(r) is another arbitrary function of r. A compari- 
son with dust models interprets b(r) as the velocity func- 
tion for the shells. (It is important to note here that we 
could have also chosen v{t 7 r) of the form v' — a(r, v) iV R' , 
where a(r,v) is again an arbitrary and suitably differen- 
ciable function, which too would allow us to integrate 
equation (0J). From equation (JHJ) let us generalize the 
form of A(t, R) as, A(t, R) = r 2 gi(r,v), where gi(r 7 v) 
is a suitably differentiable function and g\{r, 1) = g(r). 
Then we similarly have A — J A^dt = r 2 g2(r, v) and at 
the initial epoch ^(r, 1) = <?(?*)• Using cauation lllll) in 
equation©, we get, 

Rv' 

2 Pg = RA >R (p + Pr ) + 2 Pr + -^ (13) 

In general, both the density and radial pressure blow 
up at the singularity, so the above equation implies that 
the tangential pressure also blows up at the singularity. 
Writing 



b(r) = l + r 2 b (r) 



(14) 



and using equations lj9"|). ljTTl) and ljT2|) in equation©, we 

get, 



VRR = -e r2gi( - r ' v ^(l + r 2 b )Re r2 32(r,v) -R + r 3 M 

(15) 

Since we are considering collapse, we have R < 0. Defin- 
ing a function h(r, v) as, 

„2r 2 g 2 (r,v) _ i 

h(r,v) = - 2 = 2g 2 (r,v) +0{r 2 v 2 ) (16) 

and using equation^SJl in eauation (|15fl . we get after sim- 
plification, 

s/vv = -Je 2r2 ^i+g2) v b Q + e 2r2 ffi (vh(r, v) + M(r, v)) 

(17) 

Integrating the above equation we have, 
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Note that the variable r is treated as a constant in the 
above equation. Expanding t(v,r) around the center, we 
get, 

t(v,r) =t(v,0) +rX(v) + 0(r 2 ) (19) 
where the function X{v) is given by, 

1 f 1 y/v(biv + vhi(v) +Mi(v)) 
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dv w y 1 ly ' 1V 'i (20) 

(b 00 v + vh (v) + M (v)) 2 



with, b 00 = b o (0),M Q {v) = M(0,v),h o = h(0,v),h = 
b' (0), Mi(«) = M, r {0, v), hi = h, r (0, v). 

From the above, the time when the central singularity 
develops is given by 



y/vdv 



y/boQV + vh (v) + Mq(v) 



(21) 



The time for other shells to reach the singularity can be 
given by the expansion, 



t s (r) =t S0 +rX(0) + O(r 2 ) 



(22) 



It is now clear that the value of X(0) depends on the 
functions 60, M. and h, which in turn depend on the ini- 
tial data at t — ti and the dynamical variable v. Thus, a 
given set of initial matter distribution and the dynamical 
profiles including the velocity of shells completely deter- 
mine the tangent at the center to the singularity curve. 
Further, from equations (17-19), we get, 



y/wJ = X(v)y/b QO v + vh (v) + M (v) + 0{r 2 ) (23) 

It is now possible to see how the final fate of col- 
lapse is determined in terms of either a black hole or 
a naked singularity. This is determined by the causal be- 
havior of the apparent horizon, which is the boundary of 
trapped surfaces forming due to collapse, and is given by 
R = F. If the neighborhood of the center gets trapped 
prior to the epoch of singularity, then it is covered and a 
black hole results, otherwise it could be naked when non- 
spacelike future directed trajectories escape from it. The 
key point is to determine if there are any future-directed 
non-spacelike paths emerging from the singularity. To 
see this, and to examine the nature of the central singu- 
larity at R = 0, r = 0, consider the equation for outgoing 
radial null geodesies, 



dt 
dr 



(24) 



We want to examine if there would be any families of fu- 
ture directed null geodesies coming out of the singularity, 
thus causing a NS phase. 

The singularity occurs at v(t s (r),r) = 0, i.e. 
R(t s (r),r) = 0. Therefore, if there are any future di- 
rected null geodesies, terminating in the past at the sin- 
gularity, we must have R — > as t — > t s along these 
curves. Now writing the null geodesies equation in terms 
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of the variables (u — r a ,R), choosing a = |, and using 
equation© we get, 



(25) 



If the null geodesies terminate at the singularity in 
the past with a definite tangent, then at the singular- 
ity we have 4^ > 0, in the (u,R) plane with a finite 
value. Hence all points r > on the singularity curve are 
covered since F/R -> oo with H -> -oo. The central 
singularity could however be naked. Define the tangent 
to the outgoing null geodesic from the singularity as, 



,. ,. R dR 
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du 



Using equation and l|23|l we then get, 
xl /2 = (3/2)V^(0)*(0) 



(26) 



(27) 



and the equation of null geodesic emerging from the sin- 
gularity is R = xqu, or in (t, r) coordinates it is given by 



t-t 3 (0) =x r 5 / 3 



(28) 



It follows now that if X(Q) > 0, then xq > 0, and we get 
radially outgoing null geodesies coming out from the sin- 
gularity, giving rise to a naked central singularity. How- 
ever, if X(0) < we have a black hole solution, as there 
will be no such trajectories coming out. If X(Q) — then 
we will have to take into account the next higher order 
non-zero term in the singularity curve equation, and a 
similar analysis has to be carried out by choosing a dif- 
ferent value of a. 

We note that in the above the functions h and M. are 
expanded with respect to r around r = and the first- 
order terms are considered. At times, however, these 
are assumed to be expandable with respect to r 2 , and 
it is argued that such smooth functions would be phys- 
ically more relevant. Such an assumption comes from 
the analyticity with respect to the local Minkowskian co- 
ordinates (see e.g. [5]), and it is really the freedom of 
definition mathematically. We may remark that the for- 
malism as discussed here would work for such smooth 
functions also, which is a special case of above. 

We thus see how the initial data in terms of the free 
functions available determine the BH/NS phases as col- 
lapse end states, because ^(0) is determined by these ini- 
tial and dynamical profiles as given by eauation (|2(Jfl . It is 
clear, therefore, that given any regular density and pres- 
sure profiles for the matter cloud from which the collapse 
develops, we can always choose velocity profiles so that 
the end state of the collapse would be either a naked sin- 
gularity or a black hole, and vice- versa. It is interesting 
to note here that physical agencies such as the spacetime 
shear within a dynamically collapsing cloud could nat- 
urally give rise to such phases in gravitational collapse 



[6]. In other words, such physical factors can naturally 
delay the formation of apparent horizon and the trapped 
surfaces. 

We believe numerical work on collapse models may 
provide further insights into these interesting dynam- 
ical phenomena, and especially when collapse is non- 
spherical, which remains a major problem to be consid- 
ered. Much numerical work has been done in recent years 
on spherical scalar field collapse [7] , and also on some per- 
fect fluid models [8] . While we have worked out explicitly 
here the emergence of null geodesies from the singularity, 
thus showing it to be naked (or otherwise) in an analytic 
manner, the numerical simulations generally discuss the 
formation or otherwise of trapped surfaces and apparent 
horizon, and these may also break down closer to the 
epoch of actual singularity formation. In that case this 
may not allow for actual detection of BH/NS end states, 
whereas important insights on critical phenomena and 
dispersal have already been gained through numerical 
methods. Probably, a detailed numerical investigation 
of the null geodesies in collapse may provide further in- 
sights here. 

As stated above, we work here with type /matter fields, 
which is a rather general form of matter which includes 
practically all known physically reasonable fields such as 
dust, prefect fluids, massless scalar fields and so on. How- 
ever, it is important to note that suitable care must be 
taken in interpreting our results. While we have shown 
that the initial data and dynamical evolutions chosen do 
determine the BH/NS end states for collapse, the point 
is, actually, all these dynamical variables are not explic- 
itly determined by the initial data given at the initial 
epoch (note that v plays the role of a time coordinate 
here). Hence these functions are determined only as a 
result of time development of the system from the ini- 
tial data if we have the relation between the density and 
pressures, that is an 'equation of state'. In principle, it 
is possible to choose these functions freely, only subject 
to an energy condition, and after that one may calculate 
the energy density, and the radial and tangential pres- 
sures for the matter. However, in that case the resultant 
'equation of state' could be quite strange in general. Of 
course, presently we have practically very little idea on 
what kind of an equation of state should the matter fol- 
low, especially closer to the collapse end states, where 
we are already dealing with ultra-high energies and pres- 
sures. Hence if we allow for the possibility that we could 
freely choose the property of the matter fields as above, 
or the equation of state, then our analysis is certainly 
valid. In such a case, however, the chosen 'equation of 
state' will be in general such that the pressures may ex- 
plicitly depend not only on the energy density, but also 
on the time coordinate. 

All the same, it is worth pointing out that the analysis 
as given above in fact does include several well-known 
classes of collapse models and equation of states as we 
show below. It has already been stated that we have a 
choice of two free functions available here, among which 
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we are choosing only one, namely, the form of metric 
function i/(t,r). Therefore, in addition to the Einstein 
equations, if any equation of state of the form p r = f(p) 
(or alternatively, pg = g(p)) is given, then it is clear 
from equation J5J that there would be a constraint on 
the otherwise arbitrary function M, specifying the re- 
quired class, if the solution of the constraint equation 
exists. Then the value of pg (or p r ) in terms of p(t,r) is 
determined from equation ©, and the analysis for the 
BH/NS phases as above still goes through. For example, 
for dust collapse models we have 

Pr =pg = (29) 

and the constraint equation gives, 

M(r,v) = M(r) (30) 

Next, for the collapse with a constant (or zero) radial 
pressure, but with a variable tangential pressures al- 
lowed, we have, from the constraint equation, 

M(r,v) = f(r) - kv 3 (31) 

where, k is the value of the constant radial pressure. The 
tangential pressure will then be given by , 

pg = k + ^A. R R(p + k) (32) 

Along with these two well-known models, the above 
analysis works for any other models also in which one 
of the pressures is specified by an equation of state, and 
which permits a solution to the constraint equation on 
M. Hence it follows that the considerations above pro- 
vide an interesting framework for the study of dynamical 
collapse, which is one of the most important open prob- 
lems in gravity physics today. 

We note here that in the consideration above, we have 
dealt with the local visibility or otherwise of the sin- 
gularity, which is the question of strong cosmic censor- 
ship. Such a locally naked singularity may or may not 



be globally visible, which could communicate to faraway 
observers in the spacetime (thus violating the weak cos- 
mic censorship) . This is actually decided by the global 
behaviour (for large values of r) of the concerned func- 
tions such as the mass function and others. Typically 
it is observed (see e.g. [1]) that once the singularity is 
locally visible, it is possible to choose the behaviour of 
functions away from the center in such a manner that 
it would be globally naked as well. A choice otherwise 
could make all the non-spacelike trajectories coming out 
of the singularity to fall back into the apparent horizon, 
thus giving only a locally naked singularity. 



The basic result we proved is, given a regular initial 
data in terms of the density and pressure profiles at the 
initial epoch from which the collapse develops, there are 
sets of dynamical evolutions, i.e. classes of solutions to 
the Einstein equations, which can evolve the given initial 
data to produce either a black hole or a naked singularity 
as the collapse end state. The exact outcome in terms 
of BH/NS end states depends on the choice of rest of 
the free functions available, such as the velocities of the 
collapsing shells. In other words, the initial data space 
can be divided into distinct subspaces, those that evolve 
into black holes, and others into a naked singularity. In 
a sense this brings out the stability of these collapse out- 
comes with respect to perturbations in the initial data, 
and the choice of matter that collapses. This would be 
useful information on the important issue of genericity 
for collapse outcomes. Of course stability theory in gen- 
eral relativity is a complex issue and what is needed is 
a study of non-spherical collapse. While some insights 
have been provided earlier by specific examples, the gen- 
eralization as given here may be useful in view of the 
significance of this problem in black hole physics, and for 
relativistic astrophysics in general. Further details shall 
be presented elsewhere. 
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